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The BH mass of nearby QSOs: a comparison of the bulge 
luminosity and virial methods 
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ABSTRACT 

We report on the analysis of the photometric and spectroscopic properties of a sample 
of 29 low redshift [z < 0.6) QSOs for which both HST WFPC2 images and ultraviolet 
HST FOS spectra are available. For each object we measure the R band absolute 
magnitude of the host galaxy, the Civ(1550 A) hne width and the 1350 A continuum 
luminosity. From these quantities we can estimate the black hole (BH) mass through 
the Mbh ~ -^buigc relation for inactive galaxies, and from the virial method based 
on the kinematics of the regions emitting the broad lines. The comparison of the 
masses derived from the two methods yields information on the geometry of the gas 
emitting regions bound to the massive BH. The cumulative distribution of the line 
widths is consistent with that produced by matter laying in planes with inclinations 
uniformly distributed between '^10° and '^50°, which corresponds to a geometrical 
factor / ^ 1.3. Our results are compared with those of the literature and discussed 
within the unified model of AGN. 

Key words: galaxies: active - galaxies: nuclei - quasars: general. 
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1 INTRODUCTION 

Nowadays there is wide consensus that most if not all nearby 
early type galaxies host a supermassive black hole (SMBH) 
at their centers (see for a recent review Ferrarese 2006) . How- 
ever a direct determination of the BH mass (Mbh) is possible 
only for a limited number of inactive galaxies in the local 
Universe, from dynamical methods (virial theorem) based on 
the observations of the kinematics of stars orbiting around 
the SMBH. An extension of this method to higher redshift 
objects is viable only for active galaxies, where Mbh can 
be derived from the dynamics of gas gravitationally bound 
to the SMBH: this is the usual case of regions emitting the 
broad emission lines (BLR). 

At low redshift correlations between Mbh and the 
global properties of the host galaxies, namely the luminosity 
(ibuigc) and the central velocity dispersion of the bulge, have 
been found. Locally therefore the measure of the properties 
of the host galaxy is an alternative method for evaluating 
Mbh. 

The aim of this paper is to compare the mass determi- 
nations through the virial method and through the global 
parameters of the hosts of close-by quasars. For QSOs it is 
very difhcult to measure the velocity dispersion of the stars 
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in the bulge, so one must rely on the Mbh — ibuigo rela- 
tionship. For nearby ellipticals, there is increasing evidence 
that this relation has the same behavior for both inactive 
and active galaxies (see for example Merritt & Ferrarese 
2001; McLure & Dunlop 2001, 2002; Bettoni et al. 2003), 
notwithstanding a substantial increase of the dispersion in 
the latter case, possibly due to the greater difficulty in the 
determination of all the parameters. Moreover, AGN hosts 
and inactive galaxies appear substantially the same objects, 
as indicated by direct morphological studies (e.g. Dunlop et 
al. 2003). This is consistent with a picture where the AGN 
phenomenon, which is linked to the gas accretion on the cen- 
tral BH, would be a phase characterizing the life of normal 
galaxies in particular during major mergers (see for exam- 
ple KaufFmann & Haehnelt 2000). For these reasons, in this 
work we start with the assumption that the correlation es- 
tablished for nearby galaxies is also valid for higher redshift 
QSOs, provided that host galaxy luminosity is corrected for 
galactic evolution to 2: ~ 0. 

A second set of QSO BH mass determinations can be 
obtained through the virial theorem, provided that an esti- 
mate of the velocity ublr and radius iiBLR of the BLR is 
available. The BLR velocity can be derived from the FWHM 
of the broad emission lines if one assumes a value for the ge- 
ometrical factor /, which is an unknown parameter of order 
unity linked to the geometry of the BLR. The measurement 



2 M. Labita, A. Treves, R. Falomo and M. Uslenghi 



of -Rblr may require cumbersome procedures, like reverber- 
ation mapping (Ulrich et al. 1984; Blandford & McKee 1982; 
Peterson 2001), that uses the continuum source variability 
to evaluate time delays between the central emission and 
the BLR gas response; an alternative is to infer the BLR 
size from the ionizing luminosity with which it is correlated 
(Kaspi et al. 2000; Vestergaard 2002; Plan, Falomo & Treves 
2005; Kaspi et al. 2005; Vestergaard & Peterson 2006). BH 
mass determinations can be then derived from a measure of 
the broad line widths and of the ionizing continuum inten- 
sity, once the value of the geometrical factor is fixed. 

In the next section we focus on the quasar sample choice 
and on the data collection and analysis ll2.1H2.2t . In section|3 
we describe the formulae used for applying the two methods 
of mass determination II3.1H3.2II and report on the results 
l|3.3|l . This analysis constrains the value of the geometrical 
factor. The issue of the BLR geometry and orientation is 
further discussed in paragraph 13.41 In section |1] our results 
are compared with those of literature and discussed within 
the unified model of AGN. 

Throughout this paper, we adopt a concordant cosmol- 
ogy with Ho = 72 km s"^ Mpc"\ fi™ = 0.3 and = 0.7. 
We have converted the results of other authors to this cos- 
mology when adopting their relations and using their data. 
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Figure 1. The redshift distribution of the 29 QSOs with both 
HST images and W spectra. 



2 THE SAMPLE AND DATA ANALYSIS 

We selected^ all QSOs from the Veron-Cetty & Veron (2003) 
catalogue, for which HST images and spectra are available. 
We further require that the images have exposure times 
longer than 1000 s and that spectra cover the range 1400- 
1700 A (JJV band). In addition we selected only objects 
hosted by ellipticals, so that the bulge component substan- 
tially coincides with the whole galaxy. 

The resulting sample includes 29 objects at z <0.6, 17 of 
which are radio loud (RLQ) and 12 are radio quiet (RQQ). 
The main properties of the objects are listed in Table 1. 
The sample is essentially uniformly distributed in redshift 
(see Fig. nj. 



2.1 HST WFPC2 images 

For 27 sources in the sample a measurement of the host 
galaxy magnitude and morphology is available in the litera- 
ture; while for 2 objects the determination of the parameters 
of the host galaxy is obtained in this work: the images re- 
duction and analysis is discussed in Appendix A and the 
relevant data are reported in Table 5. In order to obtain a 
homogeneous set of absolute magnitudes of the host galax- 
ies, we collect from literature the values of apparent magni- 
tude (see Table l^l for references). Since Floyd et al. (2004) 
didn't report the values of apparent magnitude for their ob- 
jects, we derived them from absolute values and corrections 
they applied. For all the objects in our sample we apply 
color transformations to the K band (following Fukugita, 
Shimasaku & Ichikawa 1995), Galactic extinction (following 
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Schlegel, Finkbeiner & Davis 1998^^) and fc-correction (fol- 
lowing Poggianti 1997) , assuming that the host galaxies are 
dominated by old stellar population (see Table|5J. If we sup- 
pose that the host gala^xies have a significant fraction of a 
young (< 1 Gyr) population, the color transformation and 
fc-correction would produce magnitudes that are systemati- 
cally smaller by ~ 0.15 mag. Anyway, there is not strong evi- 
dence that QSO hosts in the nearby Universe {z < 0.5) have 
a significative young stellar component (Nolan et al. 2001). 

In Table 121 we report the i?-band absolute magnitudes 
(Mr) corrected for passive evolution (following Bressan, 
Chiosi & Fagotto 1994) for the 27 objects discussed above 
and for the 2 measured by us. Note that we have consid- 
ered only QSOs hosted by elliptical galaxies; however for 
2 objects (1001-1-291 and 1444-1-407) there is disagreement 
on the morphology of the host galaxy in literature, and our 
choice of an elliptical morphology follows respectively the 
arguments of Floyd et al. (2004) and BahcaU et al. (1997). 
For 19 objects, more than one estimate of the apparent mag- 
nitude was available. In these cases we have considered the 
mean value. The average values of the absolute R magni- 
tudes (not corrected for passive evolution) are: 

(Mr) = -23.47 ± 0.81 QSOs; 
(Mr) = -23.68 ± 0.44 RLQs; 
(Mfl) = -23.18 ± 1.10 RQQs 

where the uncertainty is the standard deviation. The gap 
between the RLQ and RQQ host magnitudes is consistent 
with the results of other researchers (see for example Hamil- 
ton et al. 2002). 



^ This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Lab- 
oratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. 
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Table 1. QSOs at z < 0.6 with HST images and spectra. These data are taken from the Veron— Getty & Veron catalogue. 
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Name 


Alternative Name 


R.A. 


(J2000.0) 


Decl (J2000.0) 


z 


Rt-o 


V 


(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


(8) 


0100+020 


PHL 0959 


01 


03 


12.9 


+02 21 10.3 


0.393 


Q 


16.39 


0133+207 


3C 047 


01 


36 


22.4 


+20 57 14.7 


0.425 


L 


18.10 


0137+012 


PHL 1093 


01 


39 


57.3 


+01 31 46.2 


0.272 


L 


17.07 


0624+691 


HS 0624+6907 


06 


30 


08.6 


+69 05 40.0 


0.370 


Q 


14.20 


0850+440 


US 1867 


08 


53 


34.2 


+43 49 01.0 


0.513 


Q 


16.40 


0903+169 


3C 215 


09 


06 


31.9 


+16 46 11.5 


0.412 


L 


18.27 


1001+291 


TON 28 


10 


04 


06.1 


+28 55 19.2 


0.330 


Q 


15.50 


1136-135 


PKS 1136-13 


11 


39 


10.7 


-13 50 43.6 


0.554 


L 


16.17 


1150+497 


4C +49.22 


11 


53 


22.3 


+49 30 21.4 


0.334 


L 


17.10 


1202+281 


PC 1202+281 


12 


04 


42.1 


+27 54 11.0 


0.165 


Q 


15.60 


1216+069 


PC 1216+069 


12 


19 


23.1 


+06 38 26.8 


0.331 


Q 


15.65 


1219+755 


MRK 205 


12 


21 


45.9 


+75 19 06.5 


0.071 


Q 


15.24 


1226+023 


3C 273 


12 


29 


06.7 


+02 03 08.6 


0.158 


L 


12.90 


1230+097 


LBQS 1230+0947 


12 


33 


28.8 


+09 31 04.9 


0.415 


Q 


16.15 


1302-102 


PC 1302-102 


13 


05 


33.0 


-10 33 19.4 


0.286 


L 


15.20 


1307+085 


PC 1307+085 


13 


09 


47.0 


+08 19 48.9 


0.155 


Q 


15.10 


1309+355 


TON 1565 


13 


12 


16.3 


+35 14 36.7 


0.184 


L 


15.64 


1416-129 


PG 1416-129 


14 


19 


05.7 


-13 10 56.5 


0.129 


Q 


16.10 


1425+267 


TON 202 


14 


27 


33.6 


+26 32 52.9 


0.366 


L 


15.68 


1444+407 


PG 1444+407 


14 


46 


45.9 


+40 35 05.0 


0.267 


Q 


15.70 


1512+370 


PG 1512+370 


15 


14 


39.2 


+36 50 37.7 


0.370 


L 


16.27 


1545+210 


3G 323.1 


15 


47 


43.5 


+20 52 16.7 


0.264 


L 


16.70 


1641+399 


3G 345 


16 


42 


58.8 


+39 48 37.0 


0.594 


L 


15.96 


1704+608 


3G 351 


17 


04 


41.4 


+60 44 30.5 


0.371 


L 


15.28 


1821+643 


HB89 1821+643 


18 


22 


02.8 


+64 20 05.3 


0.297 


Q 


14.10 


2128-123 


PKS 2128-12 


21 


31 


35.3 


-12 07 04.8 


0.501 


L 


16.11 


2135-147 


PKS 2135-14 


21 


37 


45.2 


-14 32 55.8 


0.200 


L 


15.50 


2141+175 


HB89 2141+175 


21 


43 


35.5 


+17 43 48.0 


0.213 


L 


15.73 


2201+315 


4C +31.63 


22 


03 


15.0 


+31 45 38.3 


0.295 


L 


15.58 



2.2 HST FOS spectra 

In order to have a uniform set of the spectral parameters, 
we performed systematically the following procedures to the 
available data. When more than one spectrum was avail- 
able for a single object, different observations were combined 
into a single spectrum. Then, the continuum was fitted with 
a power low, removing the regions contaminated by spec- 
tral lines and by the Fell or Hen emissions. The rest-frame 
Civ lines of the contirmum subtracted spectra are shown in 
Fig-H 

In order to evaluate the virial velocity of the broad line 
region from UV spectra, we measured Civ line-widths, fol- 
lowing i.e. Plan et al. 2005. The line widths have been mea- 
sured by fitting Civ emission lines, excluding the line regions 
contaminated by absorption features or by the Hell(1640A) 
emission. The lines have been fitted with the sum of two 
Gaussians with the same central wavelength, and then we 
derived the FWHM of the sum function. This is just a math- 
ematical tool, and it is not related to a physical decompo- 
sition into components of different origin. We tested that, 
inserting a narrow (< 1500 km/s) Civ component, the re- 
sults remain substantially unchanged, consistently with the 
numerous evidences that in the UV spectra of QSOs the nar- 
row line components are weak and rare (about 2 per cent 
occurrence rate, see Wills et al. 1993; Baskin & Laor 2005). 
Fig. |5| shows the adopted fits. 

Since in the literature virial BH mass estimates are of- 
ten based on H/3 measurements, it is interesting to compare 



our FWHMs of Civ with those of H/?. We found H/? mea- 
surements for 25 objects of our sample in the dataset pre- 
sented by Marziani et al. 2003: the average FWHM of H/3 
(5300 ± 2200 km s~^) is roughly consistent with that derived 
from our Civ measurements (4200 ± 1100 kms~^). Note the 
larger dispersion of the H/3 measurements that is likely due 
to the non homogeneity of the observations from which these 
data were derived. The consistency between Civ and H/3 
FWHMs is also supported by comparisons based on larger 
samples of objects (see Vestergaard et al. 2006 and refer- 
ences therein). 

To evaluate the size of BLR through Eq. |S] we need to 
provide a measure of the ionizing luminosity at 1350 A. To 
this aim, when a direct measure was not possible (3 cases), 
the value of the continuum flux (Fx) at the requested wave- 
length was extrapolated from the best-fitting power law. 
Continuum intensities were corrected for Galactic extinc- 
tion, consistently with the previous section, and luminosity 
at 1350 A was calculated for each object. 

The relevant quantities are reported in Table |3 The 
uncertainty of the continuum level depends on the S/N ra- 
tio of the observations and on the intrinsic variability of the 
sources, giving an error of ~ 10 per cent on Fx. On the other 
hand, the uncertainty on the measured FWHM is dominated 
by the possible presence of absorption features and contam- 
inating weak Hen emission nearby the Civ line. According 
to the shape of the Civ line, the errors on the FWHM range 



4 M. Labita, A. Treves, R. Falomo and M. Uslenghi 




: 0137 + 012 11 




llLlA^,wlilAllkW7 , , 





1400 1500 1600 1700 



400 1500 1600 1700 



1400 1500 


160D 1700 


1 1 1 1 1 1 1 1 1 ir 

0903+169 jl 
/ 1 


1 1 1 1 1 1 1 1 1 1' 

■ 


1400 1500 


1600 1700 


. 1150+497 1 

iUii>Ai<.ii...i,i,. .^.^J , 




1400 1500 


1600 1700 


1 1 1 1 1 1 1 1 1 II 
■ 1319+755 1 

11«A,...I..JlAi/,l 


1 1 1 1 1 1 1 1 1 1 


1400 1500 


1600 1700 


■ 1302-ioa II 

1400 1500 


1 1 1 1 1 1 1 1 1 1. 



Figure 2. HST WFPC2 rest-frame continuum-subtracted spectra: Civ lines. In the I'-axis, A is expressed in A; in the y-axis specific 
flux is expressed in units of 10~^^ erg cm^^ s^^ A^^. The best— fitting function is also plotted for each line. 



from 4 to 12 per cent, which illustrates the good quality of 
the data. 



3 BH MASS DETERMINATION 

3.1 BH mass from the A/bh — ibuige 

In nearby inactive galaxies, the correlation between Mbh 
and the luminosity of the bulge component of the host has 
been deeply investigated (e.g. Magorrian et al. 1998; Ko- 
rmendy & Gebhardt 2001; Bettoni et al 2003). Here we 
use the relationship obtained by Bettoni et al. (2003), who 
adopted our same definition and corrections for the magni- 
tudes. The relation corrected for the chosen cosmology is: 



log(AfBH/M(: 



-0.50A/ii - 2.60 



(1) 



where Ma is the absolute magnitude of the bulge component 
of the host galaxy. The rms of this fit is 0.39. This relation 
is close to that obtained by McLure & Dunlop (2002) for 
AGN host galaxies. 



Thus, from the absolute magnitude of the host galaxies, 
corrected for passive evolution (see Table |5J, we can deduce 
the BH masses using Eq. [H These masses are reported in 
column (2) of Table 0] The mean values are: 



(log(A/BH/M0)) 
(log(AfBH/M0)) 
{log(AfBH/M0)) 



8.99 ±0.38 
9.08 ±0.21 
8.86 ±0.52 



QSOs; 
RLQs; 
RQQs 



where the uncertainty is the standard deviation. 



3.2 BH mass from virial method 

Recent studies (Gaskell 1998; Peterson & Wandel 2000; Geb- 
hardt et al. 2000; Ferrarese et al. 2001) suggest that the 
black hole gravity dominates radiation pressure effects on 
the BLR. Assuming this and supposing that the system is 
in equilibrium, it is possible to use the virial theorem to 
evaluate the central BH mass: 



AfB 



G ^ -Rblr (j'blr) 



(2) 
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Figure l^l — continued. 



where 7?blr is the radius of the broad line region and vblk 
is the gas velocity at -Rblr- 

Under the hypothesis that the spectral line-width is 
caused by Doppler effect, the BLR velocity can be expressed 
as: 



l^BLR = / ■ FWHM 



(3) 



where FWHM refers to the broad emission lines produced 
in the BLR and is expressed in velocity units, and / is a 
factor related to the BLR geometry. For an isotropic field of 
velocities, / = V3/2. The issue of the value of / is discussed 
later on. 

A relation which links 7?blr to the ionizing luminosity 
at 5100 A was calibrated on reverberation measurements 
(Wandel, Peterson & Malkan 1999, Kaspi et al. 2000). Re- 
cently, this phenomenological formula has been extended 
outside the optical spectral range, giving expressions to de- 
rive -Rblr from the continuum luminosity in the UV and X- 
band (Vestergaard 2002; Pian et al. 2005; Kaspi et al. 2005; 
Vestergaard & Peterson 2006). To infer an estimate of the 
BLR radius, we adopt the _RBLR-iA(1350A) relation given 



in Pian et al. 2005 that was calibrated directly on reverber- 
ation measures of i?BLR: 



i?BLR = (22.4 ± O.i 



ALa(1350A) 



1044 erg s-i 



It — days, (4) 



which is consistent with that proposed by Kaspi et 
al. (2005), based on the study of 32 reverberation mapped 
AGN. 

Using Eq.|3and0]in the virial theorem (Eq.|5J, one can 

derive the BH mass: 



MBH = ,f^-4.4- lO'^ 



ALa(1350A) 



10-44 erg s-i 



FWHM ^ ^ 



1000 km s-\ 



Mq (5) 



Note that AIbh depends strongly on the value of the geo- 
metrical factor /. 

We used the values of ALa(1350A) and FWHM reported 
in Table |3] The corresponding mass values are calculated 
with Eq. |K| and are given in Table |11 under the assumption 
of isotropy of the BLR geometry (/ = \/3/2). 
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Table 2. Host galaxies absolute magnitudes of the 29 QSOs. 
Notes - (2) 1. Hooper, Impey & Foltz 1997; 2. Hamilton, Caser- 
tano & Turnshek 2002; 3. Floyd et al. 2004; 4. Pagani, Falomo 
& Treves 2003; 5. Dunlop et al. 2003; 6. Boyce et al. 1998; 7. 
Bahcall et al. 1997; 8. Kirhakos et al. 1999; 9. this work. - (3) 
ij-band mean apparent magnitudes of the host galaxies corrected 
for Galactic extinction. — (4) fc— correction values following Pog- 
gianti (1997). — (5) Elliptical galaxy passive evolution correction 
[mag] following Bressan et al. (1994). - (6) Host galaxies absolute 
magnitudes corrected for passive evolution. 



Name 


Authors 


Rext 


k(R) 


E{R) 


Mn 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


01 00-1-090 


1 2 


19.83 


0.47 


-0.36 


-21.90 




2 4 


19.23 


0.54 


-0.39 


-22.73 




2 5 6 


17.42 


0.28 


-0.23 


-23.34 


0fi24-l-fi91 


3 


17.01 


0.43 


-0.33 


-24.57 


08^0—1 44 


g 


19.13 


0.75 


-0.48 


-23.44 


OQO'i-Ul fiQ 


9 4 


19.07 


0.51 


-0.38 


-22.81 


1 001 -I-2Q1 


3 


17.82 


0.36 


-0.29 


-23.45 


1136-135 


9 


19.12 


0.87 


-0.52 


-23.73 


1150+497 


3 


17.59 


0.37 


-0.30 


-23.69 


1202+281 


2, 7 


17.14 


0.17 


-0.17 


-22.37 


1216+069 


2, 6 


18.96 


0.37 


-0.29 


-22.32 


1219+755 


2 


15.14 


0.06 


-0.08 


-22.57 


1226+023 


2, 7 


15.43 


0.16 


-0.16 


-23.93 


1230+097 


3 


17.96 


0.52 


-0.38 


-23.93 


1302-102 


2, 6, 7 


17.42 


0.31 


-0.26 


-23.44 


1307+085 


2, 7 


17.50 


0.15 


-0.16 


-21.85 


1309+355 


2, 4, 7 


16.51 


0.19 


-0.18 


-23.27 


1416-129 


2 


17.42 


0.12 


-0.14 


-21.44 


1425+267 


2, 8 


18.44 


0.42 


-0.33 


-23.08 


1444+407 


7 


18.01 


0.29 


-0.24 


-22.67 


1512+370 


2, 4 


18.45 


0.43 


-0.33 


-23.13 


1545+210 


2, 7 


17.60 


0.28 


-0.24 


-23.07 


1641+399 


8 


19.31 


1.01 


-0.56 


-23.82 


1704+608 


2, 6 


18.00 


0.43 


-0.33 


-23.58 


1821+643 


3 


16.50 


0.32 


-0.27 


-24.45 


2128-123 


6 


19.35 


0.71 


-0.47 


-23.14 


2135-147 


2, 5, 7 


16.94 


0.21 


-0.19 


-23.08 


2141+175 


2, 5 


17.03 


0.23 


-0.20 


-23.12 


2201+315 


2 


16.65 


0.32 


-0.26 


-24.30 



3.3 Comparison of BH masses 

In Fig. Elwe compare the BH masses obtained with the pho- 
tometric and virial methods. As clearly shown in panel (a), 
the simple case of isotropic distribution of the emitting re- 
gion yields inconsistent results: a systematic difference be- 
tween the mass determinations obtained with the two meth- 
ods is apparent ( { A log(MBH/MQ) ) = 0.36). 

Nevertheless, a linear relation fits well the data point 
and the slope of the best fit is consistent with unity (a — 
1.1 ±0.2). This suggests that a way of matching the two BH 
mass estimates is to consider the geometrical factor / as a 
free parameter, and to search for the value which minimizes 
the sum of the differences between the masses in logarith- 
mic units from the two methods. This yields / = 1.3 ± 0.1. 
Assuming this value, we report in Table|l]the corresponding 
virial BH masses. 

In panel (b) of Fig.|H]we show the comparison of the BH 
masses obtained with the photometric and virial methods, 
adopting / = 1.3. 



Table 3. Rest— frame spectroscopic data of the 29 analyzed ob- 
jects. Notes - (2) Rest-frame Civ FWHM values - (3) Rest- 
frame flux values — (4) Galactic Extinction values following 
Schlegel et al. 1998 - (5) Luminosity at 1350 A. 



Name 


FWHM 


F;^(1350A) 




log(AL;^) 






L cm~ s A J 


[mag] 




(1) 


(2) 


(3) 


(4) 


(5) 


0100+020 


3.97 


1.69 


0.17 


45.15 


0133+207 


4.48 


1.16 


0.52 


45.21 


0137+012 


6.78 


0.81 


0.23 


44.49 


0624+691 


3.58 


16.72 


0.80 


46.35 


0850+440 


4.22 


3.00 


0.32 


45.73 


0903+169 


4.99 


0.69 


0.34 


44.88 


1001+291 


3.07 


11.33 


0.18 


45.81 


1136-135 


2.94 


4.19 


0.40 


45.99 


1150+497 


3.97 


3.20 


0.17 


45.26 


1202+281 


2.56 


2.37 


0.17 


44.45 


1216+069 


2.56 


11.21 


0.18 


45.81 


1219+755 


2.94 


19.71 


0.37 


44.74 


1226+023 


3.20 


108.12 


0.17 


46.05 


1230+097 


3.33 


5.57 


0.17 


45.72 


1302-102 


3.97 


31.22 


0.34 


46.16 


1307+085 


3.33 


9.20 


0.28 


45.03 


1309+355 


3.20 


12.13 


0.10 


45.24 


1416-129 


3.84 


3.16 


0.79 


44.57 


1425+267 


8.06 


4.38 


0.15 


45.48 


1444+407 


5.12 


15.91 


0.11 


45.71 


1512+370 


3.97 


8.17 


0.18 


45.79 


1545+210 


4.22 


5.66 


0.34 


45.35 


1641+399 


4.35 


3.93 


0.14 


45.93 


1704+608 


5.12 


10.92 


0.19 


45.91 


1821+643 


3.97 


37.10 


0.34 


46.27 


2128-123 


4.73 


7.59 


0.58 


46.22 


2135-147 


4.86 


5.22 


0.42 


45.09 


2141+175 


4.22 


4.54 


0.90 


45.27 


2201+315 


4.73 


14.90 


0.99 


46.13 



3.4 Constraining the BLR geometry 

The results of the previous paragraph strongly suggest that 
the simplifying hypothesis of complete isotropy of the BLR, 
which leads to f — \/3/2, is to be revised. Therefore we 
consider a scenario where there is a component of the BLR 
velocity field that lies on the plane of the accretion disk (see 
also McLure & Dunlop 2002; Jarvis & McLure 2006), in 
order to explore the connection between the BLR geometry 
and the value of the geometrical factor. 

Within the unification model, we suppose that AGN 
are randomly oriented in space and that they are observed 
with the characteristics of QSO if they lie with an inclination 
angle 6 in a range Omin — Omax ■ The probability distribution 
to observe an object in this range is: 



Po{0) = 



(6) 



cos 6'min — cos 6'max 

and Pg{0) — outside. 

The observed FWHM of the broad emission lines is 
given by: 



FWHM = 2 ( ci--i=sin6l + C2--^ ) vblr 



(7) 



where sin^ accounts for the projection of ijblr on the line 
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Table 4. BH mass determinations. Notes - (2) Mbh from the 
bulge luminosity (Eq. 0. - (3) Mbh through the virial method 
(Eq.|3 adopting / = \/3/2. - (3) A/bh through the virial method 
(Eq.EJ adopting / = 1.3. 



Name 


log(MBH/MQ) 


log(MBH/M0) 


log{MBH/MQ) 




from Lbuige 


/ = V3/2 


/ = 1.3 


(1) 


(2) 


(3) 


(4) 


ni nn_i_non 


o.oO 


a /1 1 

0.4X 


fi TV 
O. 1 1 


ni QQ 1 on'7 


0.(0 


S 

tt.Ou 


Q OI 

o.yi 




n r\'7 

y.u ( 




Q QQ 
O.OO 


nfiO/1 _i_(^oi 


y.oy 


Q nc; 
y.uo 


O /1 1 

y .41 


na c;n_i_/i yi n 


O 1 o 

y. iz 


o.oZ 


O 1 Q 

y . lo 


nOAQ-Ll fiO 


fi fin 
o.oU 


S AA 
cS.44i 


Q fin 

O.OU 


1 nni 1 oni 


O 1 O 

y. iz 




o.yo 


1 1 1 'i^ 


y.zD 


O.DO 


o no 
y .uz 


1 1 f^n_i_/io'7 


y.zo 


a AQ 

O.40 


o.o4 




o.Oo 


( -DU 


( .yo 


1 oi 1 nfio 


o.Oo 


Q /I Q 

cS.4o 


Q vn 

o. /y 


1 OI n 1 vF^c: 


o.Do 


? on 
i .yu 


o.ZD 




y.o ( 


Si 77 


O 1 Q 


1230+097 


9.36 


8.60 


8.96 


1302-102 


9.12 


9.03 


9.39 


1307+085 


8.32 


8.18 


8.54 


1309+355 


9.03 


8.28 


8.64 


1416-129 


8.12 


8.03 


8.39 


1425+267 


8.94 


9.23 


9.59 


1444+407 


8.74 


8.97 


9.33 


1512+370 


8.97 


8.80 


9.16 


1545+210 


8.94 


8.59 


8.95 


1641+399 


9.31 


8.97 


9.33 


1704+608 


9.19 


9.09 


9.45 


1821+643 


9.63 


9.09 


9.45 


2128-123 


8.97 


9.22 


9.58 


2135-147 


8.94 


8.55 


8.91 


2141+175 


8.96 


8.54 


8.90 


2201+315 


9.55 


9.16 


9.52 



of sight and the coefScients \/2 and \/3 refer to an isotropic 
distribution on the plane of the accretion disk and in the 3D 
space respectively. The parameters ci and C2 represent the 
relative fraction of FWHM due to the two velocity compo- 
nents. Note that, according to Eq. |21 the mean value of / 
can be immediately evaluated from Eq. 



where (sin ff) is averaged over the probability distribution of 
the inclination angle (Eq. EJ. With these premises, one can 
infer the expected cumulative FWHM distribution through 
Eq. © and 0. 

In Fig. m we compare the distribution of the observed 
FWHM of the Civ lines with the expected distribution, 
in the cases of the two extreme assumption of complete 
isotropy (ci = 0) and disc-like BLR geometry (c2 = 0). 

In the isotropic case we assume that the «blr has an 
intrinsic Gaussian distribution with dispersion a. The ob- 
served data are consistent with the expected distribution 
if ^p-y^n^^ = 23.4 per cent, giving a Kolmogorov-Smirnov 
compatibility test with P(KS)~ 99.5 per cent (see panel (a) 
of Fig.^J. However, this isotropic scenario implies / ~ \/3/2 
(see Eq. |H|with ci = 0) , and therefore it is inconsistent with 




log(MBH/Me) from Lj,, 



Figure 3. Virial BH mass determinations, adopting / = V^/2 
in panel (a) and / = 1.3 in panel (b), are plotted against mass 
values from the Mbh — ^bulgo relation. Filled (open) symbols 
indicate radio loud (quiet) objects. Squares highlight the objects 
with strong absorption features contaminating the Civ line. The 
solid line is the bisector; the thick dashed line is the best fit to 
the data. 



the comparison of the BH masses through the two methods 
(paragraph I3.2|l . 

In the case of disc-like geometry, we suppose that the 
FWHM dispersion is dominated by orientation (i.e. the 
intrinsic wblr dispersion is negligible). Fig. shows the 
Kolmogorov-Smirnov probability as a function of Smin and 
^MAx- Superimposed, we also plot the lines corresponding to 
fixed values of /. It is noticeable that taking Qyw ~ 5° — 20° 
and Smax ~ 45° — 50° one has P(KS)> 99.5 per cent (see 
also panel (b) of Fig. |1J and full consistency with the ge- 
ometrical factor value (/ ~ 1.3) deduced in the previous 
paragraph. 
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Figure 4. The cumulative distribution of the observed Civ 
FWHM (sohd line) is compared in panel (a) with the Gaussian 
FWHM distribution which is expected under the assumption of 
isotropic BLR geometry (dotted line). In panel (b), the same dis- 
tribution of observed FWHM (solid line) is compared with the 
expectations in the case of disc— like BLR geometry (dotted line), 
provided that inclination angles are distributed between 12° and 
48° (see text for details). 




P>50% 



Figure 5. Kolmogorov— Smirnov probability values for the com- 
patibility between the cumulative distribution of the observed Civ 
FWHM in our sample and the expected distribution constructed 
under the hypothesis of disc-like BLR geometry, in function of 
the minimum and maximum inclination angles. The darkest area 
corresponds to the highest compatibility, with P(KS)> 99.5 per 
cent. The dashed lines indicate the locus of angles pairs which 
give the same geometrical factor (see text for discussion). 



they are uniformly oriented between 5° — 20° and 45° — 50°. 
This result is in good agreement with the unified model of 
AGN (see Barthel 1989; Antonucci & Miller 1985; Antonucci 
1993; Urry & Padovani 1995). 

One possible caveat on this result is that we used Civ 
FWHM measurements instead of H/3 ones. If we adopted the 
H/3 FWHMs (see Section [2.21 . we would obtain on average 
BH masses 0.2 dex larger. This would imply a value of the 
geometrical factor / ~ 1, which is still inconsistent with the 
isotropic picture. 

The quasar BH masses range over two decades (10* — 
10^'^ Mq). Both the photometric and the virial approaches 
indicate that RLQ masses exceed IO^^ '^Mq, while RQQs are 
distributed in the whole mass range. RLQs have on average 
BH masses greater than RQQs; the difference 5 between the 
mean values of log Afan of the RLQ and RQQ subsamples 



4 SUMMARY AND CONCLUSIONS 

We have evaluated the BH masses of a sample of 29 low red- 
shift quasars using both a photometric and a virial method. 
This comparison suggests the value of the geometrical fac- 
tor, / = 1.3. We also found that the cumulative FWHM 
distribution of Civ(1550 A) is consistent with both a disc- 
like BLR geometry and an isotropic velocity field, provided 
that in the latter case there is some distribution around the 
mean value of the velocity. However, the virial BH masses 
coincide on average with the photometric ones only in the 
disc picture, which strongly supports the hypothesis that 
the broad line region velocity field lies on the plane of the 
accretion disc. In this scenario, one can also put some con- 
straints on the inclination angles of QSOs, obtaining that 



5 = {log Mbh(RLQ) ) - ( log Mbh (RQQ) ) = 0.26. 

This is in good agreement with the results obtained by Dun- 
lop et al. (2003) and Falomo et al. (2004), who found respec- 
tively S = 0.34 for 23 QSOs at (z) = 0.2 and 5 = 0.33 for 16 
QSOs at higher redshift {{z) = 1.5). 

The general agreement between the photometric and 
the spectroscopic methods may be considered as an indirect 
indication that the Mbh — ^buigc correlation obtained locally 
can be extrapolated at higher redshift {z ~ 0.5). 

Our results on the BLR geometry agree well with that 
of the seminal paper of McLure & Dunlop (2001), who pro- 
posed that / = 1.5 and obtained that a sample of 30 QSOs 
and 15 Seyfert I galaxies is well described by a disc-like 
BLR geometry, provided that it is composed of two popula- 
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tions with inclination angles uniformly distributed within 
0° < e < 20° and 0° < 6 < 46° respectively. Our set 
of quasars is comparable to theirs, but the spectroscopy 
is of much better quality. In fact, the scatter in terms of 
black hole mass around the best-fitting relation of Fig. |3 
(cr — 0.30) is appreciably smaller than that obtained by 
McLure & Dunlop {a = 0.59). Our results are in good agree- 
ment also with the findings by Onken et al. (2004) , who com- 
pared Mbh derived by the virial method with those from 
the stellar velocity dispersion for a sample of 16 low redshift 
AGN, obtaining that the geometrical factor, following our 
definition, is / ~ 1.2. 

In conclusion, we believe that there is strong evidence 
that, regardless of the method used to derive BH masses, 
the geometrical factor results in the range 1.0 < / < 1.3. 
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APPENDIX A: HST IMAGING OF 2 QSO 

In Table 5 we report the journal of observations for the two 
objects without available measurement of the host galaxy 
from the literature. A two-dimensional analysis has been 
carried out using AIDA (Astronomical Image Decomposi- 
tion and Analysis, Uslenghi & Falomo 2006, in preparation; 
see also Kotilainen et al. 2005), a software package specifi- 
cally designed to perform two dimensional model fitting of 
QSO images, providing simultaneous decomposition into the 
nuclear and host galaxy components. 

PSF modelling is by far the most critical step of 
analysing QSO host galaxies. For HST images, the software 
package TinyTim (Krist & Hook 1997) can produce a PSF 
model that is very accurate in the inner part (up to 1-2") 
but it does not properly model the external fainter halo pro- 
duced by the scattered light (e.g. Scarpa et al. 2000). To 
account for this we therefore modelled the external part by 
adding an exponential term to the TinyTim PSF. The pa- 
rameters of these components have been evaluated by fitting 
the PSF model to a number of stars in the field. 

For each QSO we have then applied a mask to exclude 
possible contamination from spurious sources, bad pixels, 
companions and other defects affecting the image. The back- 
ground was estimated from the signal computed in a circular 
annulus centered on the object and its uncertainty obtained 
from the comparison of the values computed in annuli at 
different radii. The masked QSO images were finally fitted 
with a model including both a point source and the host 
galaxy assuming a de Vaucouleurs law convolved with the 
PSF. The results of the best fit parameters are reported in 
Table 5. A fit with a disc galaxy yielded significant poorer 
results, consistently with the high luminosity of the nucleus 
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Table 5. Journal of the unpublished HST WFPC2 observations and QSOs measured parameters. Notes — (8, 9) Nucleus and host 
galaxy apparent magnitudes in the observation filter. - (10) Host galaxy elliptical radius. - (11) Host galaxy apparent magnitude in the 
J?Cousin band. 



Name 


z 


Proposal ID 


Obs. Date 


Filter 


Tcxp 


Camera 




mhost 


r'o 


R 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


0850+440 


0.513 


5949 


07 Feb. '96 


F702W 


2400 


WF3 


16.93 


18.80 


1.03" 


19.21 


1136-135 


0.554 


6619 


10 Jun. '97 


F702W 


2100 


WF3 


16.52 


18.82 


1.02" 


19.23 



0850+440 



1136-135 




2 4 6 

Radius (orcaec) 



Figure 6. QSOs images, 2-D residuals after the PSF subtraction, and radial profiles: data points (squares), PSF (dotted lines), galaxy 
model (dashed lines), best— fitting convolution functions (solid lines). 



(see Dunlop et al. 2003). In Fig.|Hlwe show the QSOs images 
and residuals. The radial profiles of the images, of the PSF, 
of the galaxy and of the best-fitting models for the 2 objects 
are also shown. 



